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We study finite size effects in a family of systems in which a parameter controls 
interaction-range. In the long-range regime where the infinite-size free energy is universal, we 
show that the finite size effects are not universal but depend on the interaction-range. The 
finite size effects are observed through discrepancies between time-averages of macroscopic 
variables in Hamiltonian dynamics and canonical averages of ones with infinite degrees of 
freedom. For a high energy regime, the relation to a pair of the discrepancies is theoretically 
predicted and numerically confirmed. We also numerically show that the finite-size effects of 
macroscopic variables in the canonical ensemble are close to ones in the dynamical systems. 

§1. Introduction 

In many body systems, the range of interactions governs physical features. Long- 
range interactions induce a phase transition, and typical ones are the mean-field in- 
teractions. For instance, the Hamiltonian mean-field (HMF) model which consists 
of planar rotators shows a second-order phase transition between a low energy or- 
dered phase and a high energy disordered phase. ^ On the other hand, short-range 
interactions are hard to occur the phase transition, particularly in low-dimensional 
systems. If the planar rotators are on a one-dimensional lattice with nearest-neighbor 
interactions, then no phase transition occurs. 

To investigate how the range of interactions affects the existence of a phase 
transition, the HMF model is extended to a model having algebraically decaying 
interactions, called the a-XY model, described by the following Hamiltonian 

i=i i=i j=i j 

Each rotator is placed on a site of the periodic d-dimensional simple lattice, that 
is N = L d for the lattice of size L, and is the shortest distance between the 
sites i and j. The a-XY model includes the HMF model as a = 0, and the nearest 
neighbor interaction realizes in a — > oo. The factor N is introduced for recovering 
the extensivity of the system, but the factor does not make the system additive 
for < a/d < 1. In this article, N is defined as ^2iLi ^/ r ?j-^ 111 the long-range 
regime < a/d < 1, Campa et al. show that the a-XY model including the external 
magnetic field essentially has the same canonical partition function with the HMF 
model in N — > oo, 3 ) and the critical value of the energy for a second-order phase 
transition is common as U = E/N = 0.75, with E being the value of the Hamiltonian. 
Recent studies of the a-XY model are briefly summarized in Ref.4). 
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In this article, in the long-range regime, we investigate how macroscopic variables 
for the a-XY system with finite N approaches the thermodynamic limit (N — ► 
oo). They could go towards canonical averages as N increases (see Appendix, and 
Refs. 5)-7) for the HMF model), but it is not obvious whether the convergence 
speed depends on a/rf or not. Some non-universalities appear as dynamical aspects 
in < a/d < l, 4 ) and hence we compare canonical averages of the macroscopic 
variables in N — > oo and time-averaged variables with finite N. The time-averaged 
variables are obtained by integrating the canonical equation of motion yielded from 
the Hamiltonian (jl-lj) . We set d = 1 for the convenience of numerical computations, 
and show that, for the supercritical energy regime, the time-averages of macroscopic 
variables with finite N algebraically approach the canonical averages taken in iV — > 
oo as N increases, and the exponents depend on a. 

§2. Theoretical Prediction 

As macroscopic variables we observe, we adopt modulus of the magnetization 
and the temperature. The magnetization vector M := (M x , My) and its modulus M 

are defined as M := X^i( cos 9«> sm( 7«)/-^ an< ^ M := J M% + Mjj , respectively. The 

temperature in the sence of dynamics is defined as := Ylj=iPj/N- Note that we 
numerically compute averages through Hamiltonian dynamics instead of canonical 
ensemble, then a fixed variable is neither T nor Tpj, but the internal energy U = E/N. 
In this section, we estimate the following two discrepancies 

5m(N):=(M) n -(M) oq , 5 T (N):=T N -T, (2-1) 

where (X) N and (X) denote the canonical average of X with finite degrees of 
freedom and one in N — > cxd respectively. (U) N and (M) N are calculated in Ref 
8) with the aid of the saddle-point approximation. 3 ) Tjy is the temperature for the 
system with iV degrees of freedom and is defined as 

U = (U) N (T N ). (2-2) 

The relation (J2-2j) bridges between Hamiltonian dynamics and canonical statistics. 
We require that the relation ()2-2|) holds in any degrees of freedom and hence (U) N (Tjv) = 
(U)oo The average (U)^ (T) is related to T and {M}^ as 

2([/) 00 (T)=r+l-[(Af) 00 (T)] 2 . (2-3) 

The relation (|2-3f) is straightforwardly obtained for the HMF model, and holds in 
the long-range regime < a < 1 due to the universality of the canonical partition 
function. 3 ^ 

Although we do not show the derivation, we can show the following relation for 
T > T c , T c being the canonical critical temperature 8 ^ 

5 T (N) = [5 M (iV)] 2 + WV~ (1 - a) . (2-4) 

The constant b represents the diagonal elements of the matrix (1/r? ), and can be 
chosen arbitrarily. We select b = 2.0 to obtain a fully positive spectrum which 
enables us to apply the Hubburd-Stratonovitch transformation. 3 ) 
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§3. Numerical Results 



Let us check the validity of our theoretical predictions of St and 5m numerically. 
The time-averaged temperature \T^\ and the order parameter (M) t are defined as 
follows: 



-N 



1 f to 



T d N {t)df (M) t (t ) 



1 f to 

— / M(t) dt. 
to Jo 



(3-1) 



Numerical integrations of the canonical equation of motion are performed by using 
a fourth-order symplectic integrator. 9 ) Thanks to the Fast Fourier Transformation 
technique, the calculation times are of order iV In N. The time slice of the integrator 
is set at 0.2 and it suppresses the relative energy error AE/E ~ 1(T 5 . Our initial 
condition is, for all values of a, the canonical equilibrium distribution of the HMF 
model, that is p and q are randomly taken from the distributions being proportional 
to exp(— (3p 2 /2) and exp(/3 {M)^ cos q) 7 ^ respectively, where j3 = 1/T. The sum of 
momenta is an integral of motion, and we initially set it to zero. The computing 
time, to in Ea. (|3-l|) . are set to more than 10 6 to realize equilibrium states in the 
sense of dynamics, where we define such equilibrium as when the time-averaged 
macroscopic variables are almost stational. Typical time sequence of both (M) t (to) 
and (Tjy) (to) ar e shown in Fig^as functions of computing time to- in Fig^ we 
observe finite size effects 5 m and St, which are discrepancies from (M)^ and T 
respectively. 




10 z 10 J l(f 
computing time t 



10° 



Fig. 1. Time- averaged values (T$)t(to) and (M) t (to) as functions of computing time to- Parameters 
are taken at N = 256, a = 0.3,0.6,0.9 and U = 0.8, for which (M) x = 0.0. The horizontal 
solid line represents the corresponding canonical temperature T = 0.6. 

For the supercritical energy regime, we investigate the finite size effects for two 
values of energy. One is the high enough energy U = 5.0, and the other one is 
supercritical but the middle energy U = 0.8. The discrepancies 5m(N) and 5t(N) 
are shown in Fig(21as functions of N. It is shown that the discrepancies algebraically 
decrease as iV increases, that is 



S M (N) oc N~^ a ' u \ 8 T (N) oc ]\rM«.tf). 



(3-2) 
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Fig. 2. Dependence on N of the two discrepancies (a) {M) t -{M) oo and (b) (T%) t -T at U = 5.0. 
Asymptotic power-law decays are observed. 



According to the theoretical prediction 1)2-4(1 . the exponent /*m should be related to 
the other (jlt as 

fj,T = min(2//M, 1 — ot). (3 - 3) 

Let us check the relation 1(3-3(1 by observing the exponents (jlm and Ht as functions 
of a. 

For the high enough energy U = 5.0 (see FigE^a)), the theoretical curve is in 
good agreement with numerical results in the interval < a < 0.6. In addition, the 
theoretical prediction is qualitatively good for U = 0.8 (see FigEfb)). For further 
discussions, we refer the reader to Ref.8). 
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Fig. 3. Exponents \lm and \xt as functions of a. (a) U = 5.0. (b) (/ = 0.8. Dashed curves represent 
theoretical predictions. 



§4. Conclusions 

In this article the a-XY model 1(1-1(1 with d = 1 has been studied for the 
long-range regime < ot, < 1, so that we clarify the role of interaction length in 
Hamiltonian systems with many degrees of freedom. 8 ) We have studied dependences 
on degrees of freedom for time-averaged macroscopic variables calculated through 
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Hamiltonian dynamics with finite degrees of freedom, contrasted with the canonical 
averages of ones in the thermodynamic limit. A non-universal behavior in < a < 1 
is observed in the asymptotic behavior of the discrepancies for the supercritical 
regime. For the supercritical regime, the discrepancy between the two averages de- 
cays as a power-type function, and the exponents depending on a are explained by 
applying the canonical ensemble with finite degrees of freedom with the aid of the 
saddle-point approximation. This non-universal behavior sheds light on the study of 
Hamiltonian systems with long-range interactions. 
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Appendix 

In this appendix, to show that the values of the canonical averages of macro- 
scopic variables are close to the values of time- aver ages of ones, we numerically study 
equilibrium values of macroscopic variables for models with finite N and d = 1 in 
each situation. In the canonical ensemble and in the dynamical system, we use a 
conventional canonical Monte Carlo algorithm and a fourth symplectic integrator re- 
spectively. First, for the model with a = 0.3, and this a is close to a = 0, Figs|I]and 
|S]show that the relation (M) N -(U) N in the canonical ensemble are close to {M) t -U 
in the dynamical system. Furthermore, it is suggested that the time-averaged values 
of macroscopic variables go towards the canonical averages of ones evaluated in the 
thermodynamic limit as N increases. Next, for the model with a = 0.9, we show 
the values of the magnetization and the temperature as functions of the energy den- 
sity in Figs El and [7| The relation (M) N -{U) N in the canonical ensemble are close 
to (M) t -U in the dynamical system again. However, in this case, the prediction 
from the canonical average evaluated in the thermodynamic limit is not close to the 
relations (M) N -(U) N and (M) t -U. 

Thus, from the observations in the two cases a = 0.3 and a = 0.9, we conclude 
that, at a fixed value of a provided a < 1, the relations {M) N -{U) N in the canonical 
ensemble are close to the relations {M) t -U in the dynamical system. Meanwhile, at 
any fixed value of N, the discrepancies between (M) N -(U) N and (Af) 00 -(^ r ) 00 , and 
between (M) t -U and (Af) 0O -(^) oo! become larger as a increases up to unity. This 
tendency of the discrepancy is consistent with ones in Fig|3] for the supercritical 
energy regime. 
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Fig. 4. Dependence on N of the values of macroscopic variables in the canonical ensemble with 
a = 0.3, (a) (M) N -{U) N and (b) T-(U) N . 
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Fig. 6. Dependence on N of the values of macroscopic variables in the canonical ensemble with 
a = 0.9, (a) (M) N -(U) N and (b) T-(U) N . 
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Fig. 7. Dependence on N of the values of macroscopic variables in the dynamical system with 
a = 0.9, (a) (M) t -U and (b) (T^) t -U. Here the lines denote the canonical averages of M and 
Ti in N -> oo. 



